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■ ABSTRACT 

O I It has been suggested that diffuse ionized gas can extend ah the way to the end of 

' . the Hi disc, and even beyon d, such as in the case of the warped galaxy NGC 253 

2 ■ ([Bland-Hawthorn et alJ[l997[ ). Detecting ionized gas at these radii could carry sig- 

^ \ nificant implications as to the distribution of dark matter in galaxies. With the aim 

^ . of detecting this gas, we carried out a deep Ha kinematical analysis of two Sculptor 

group galaxies, NGC 247 and NGC 300. The Fabry-Perot data were taken at the 36- 
cm Marseille Telescope in La Silla, Chile, offering a large field of view. With almost 20 
hours of observations for each galaxy, very faint diffuse emission is detected. Typical 
emission measures of 0.1 cm~^ pc are reached. For NGC 247, emission extending up 
to a radius comparable with that of the Hi disc (r ~ 13') is found, but no emission is 

■ seen beyond the Hi disc. For NGC 300, we detect ionized gas on the entirety of our 
^vq \ field of view (rmax ^ 14'), and find that the bright Hll regions are embedded in a dif- 
fuse background. Using the deep data, extended optical rotation curves are obtained, 
as well as mass models. These are the most extended optical rotation curves thus far 
for these galaxies. We find no evidence suggesting that NGC 247 has a warped disc, 
and to account for our non detection of Ha emission beyond its H I disc, as opposed 
to the warped galaxy NGC 253, our results favour the model in which only through a 
warp, ionization by hot young stars in the central region of a galaxy can let photons 

■ escape and ionize the interstellar medium in the outer parts. 

^ ■ Key words: galaxies: NGC 247, NGC 300 - galaxies: kinematics and dynamics - 

- - ■ galaxies: ISM - instrumentation: interferometers - techniques: radial velocities. 
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1 INTRODUCTION 

It has usually been thought that the H i disc in spiral galaxies 
extends to radii well beyond the optical disc. Hi rotation 
curves have therefore led to extensive studies of the dark 
matter distrib ution at large radii. 

However, iBland -Hawthorn et al. 1 (|l997l ) suggested that 
diffuse optical emission could in principle be found at radii 
larger than the Hi disc. They initially suggested that be- 
yond a certain radius, cold gas could no longer support it- 
self against ionization by the ambient radiation field. To test 
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this theory, they obtained very deep optical spectra for one 
of the brightest Sculptor group galaxies, NGC 253, and suc- 
ceeded in detecting ionized gas (Ha and [Nil]) beyond the 
H I disc. The rotation curve they derived also showed signs 
of a possible decline. However, their results led them to con- 
clude that the existence of this extended ionized gas was not 
due to the ambient radiation field, but more likely due to 
hot young stars in the central regions of the galaxy ionizing 
the outer disc through a strong warp present in the galaxy's 
disc. 

Detecting diffuse ionized gas can therefore not only pro- 
vide insight as to the kinematics of a galaxy at large radii, 
but also provide strong constraints on the source of its ion- 
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ization. In this context, IPicaire et alJ (|2008l ) carried out a 
deep Ha kinematical analysis of NGC 7793, another Sculp- 
tor group galaxy. They used the 36-cm Marseille Telescope 
located in La Silla, Chile, offering a large field of view. With 
almost 20 hours of observations, they detected ionized emis- 
sion all the way to the end of the H i disc and confirmed that 
NGC 7793 had a truly declining rotation curve. In order to 
complete this study, another observing run was carried out 
in October 2007 with the objective of studying the Ha kine- 
matics and extent of the diffuse ionized gas of three other 
bright Sculptor group galaxies: NGC 253, NGC 247 and 
NGC 300 . The results concerning NGC 2 53 have been pub- 
lished in iHlavacek-Larrondo et al For this galaxy, 
we were able to find optical emission reaching out to 1.6 
times the maximum radius of the Hi disc. We succeeded 
not only in detecting ionized gas beyond the detections of 
iBland-Hawthorn et al.1 (Il997l ). but were also able to confirm 
the declining part of the rotation curve. Our observations 
also favoured an ionization source by hot young stars in the 
central regions of a warped galaxy. 

We now present the results concerning NGC 247 and 
NGC 300. For NGC 247, we seek to see if diffuse optical 
emission can also be seen on a scale similar to the Hi disc, 
and even beyond. For NGC 300, the Hi disc extends to r = 
20^ but previous Ha stu dies have only been a ble to detect 
ionized gas up to r 10^ (|Marcelin et aL 1985|l and down t o 



Table 1. Optical parameters of NGC 300 and NGC 247 



emission measures (EM) of 6 cm~^ pc (jHoopes et al.|[l996l V 
Although we will be limited by our field of view for this 
galaxy, we will be able to reach further distances (r ^ 14^) 
and at least an order of magnitude deeper in flux, allowing 
us to further study the Ha kinematics of the galaxy. 

We first present the data (Section 2), and then the ve- 
locity and total integrated Ha emission maps (Section 3). 
In Section 4, we derive the kinematical parameters and ro- 
tation curves, and we show the mass models in Section 
5. The results are discussed in Section 6, and summa- 
rized in Section 7. A distance of 1.80 Mpc for NGC 300 
and 2. 53 Mpc for NGC 247 have been adopted in this 
study (|Puche Carignanlll988l V The optical parameters of 
NGC 300 and NGC 247 are summarized in Table [H 



2 FABRY-PEROT OBSERVATIONS AND 
DATA REDUCTION 

The FP observations of NGC 253, NGC 247 and NGC 300 
were taken with the 36-cm Marseille Telescope at La Silla, 
Chile. The reduction steps are the same as those used in the 
study of NGC 253 in Hlavacek-Larro ndo et al.l (|2010l l and 
are described in this section. 

A commercial camera Andor iXon was used, equipped 
with an EMCCD (Electron Multiplying Charge Coupled De- 
vi ce) detecto r (details concerning this detector can be found 
in lDaide et' al. 2004, 2006, 2008, 2009). The Quantum Effi- 
ciency (QE) is ^ 90 per cent at Ha. The interference filter 
has a central wavelength of Ac = 6565.0 A, and a FWHM 
of 30 A. Maximum transmission was ^ 80 per cent at 20° C. 
The FP interference order was p = 765 at Ha, with a Free 
Spectral Range (FSR) of 8.58 A covering 40 channels of 
0.21 A each. We calibrated the data in wavelength with a 
neon lamp (A = 6598.95A). The EMCCD has dimensions 
of 512 X 512 pixels, and an image scale of 2.8^^ pixel" ^ (to- 
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Figure 1. Comparison between the fluxes obtained for NGC 300 
(left) and NGC 24 7 (right) in this study, and those of 
Gaustad et al ] (I2OOID . Each point represents a bright Hll region 
in a galaxy. The continuous line has a slope set to unity. 



tal field of view is 23.6^ x 23.6^). The exposure time was 
15 seconds per channel. For NGC 300, a total of 115 cycles 
were observed (1150 minutes or 19.2 hours), whereas for 
NGC 247 a total of 123 cycles were observed (1230 minutes 
or 20.5 hours). Since NGC 300 has a larger spatial coverage, 
it was also necessary to observe separate cubes to extract 
the sky background. For every 6 cycles taken for NGC 300, 
^ 2 cycles of sky were observed (by moving the telescope 
20' to the South). 

We calibrated in flux the data with the SHASSA 
(Southern Ha Sky Survey) catalogue ([Gaustad et al.|[200lL 
http://amundsen.swarthmore.edu/), but the catalogue was 
flrst corrected for [Nil] contamination (AA 6548 A, 6583 A). 
For NGC 300, the intensity of the [Nil] contamination was 
taken as 0.20 times that of Hq, and for NGC 24 7, it was 
taken as 0.24 that of Ra (jKennicutt et al.1 l2008h . Fig. □ 
compares the flux values obtained for NGC 300 (left) and 
NGC 247 (right) in this study, and those of iGaustad et al.l 
(|200lh . We estimate that our flux values are within 20 — 30 
per cent accurate, based on the s c atter between our flux 
values, and those of iGaustad et al.1 (|200ll ). 

We now outline the reduction steps. First, we corrected 
our data for variations of the eflficiency on the detector as 
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Figure 2. Field of view of NGC 247 and NGC 300 (linear grey-scale), each totalling 23.6 x 23. 6^ These are the raw data sets (summed 
over 40 channels), before applying the 5x5 pixels binning and gaussian fits to the line profiles. Regions contained within red contours 
show the affected area by the contaminating reflection, and appear less bright because of overcorrection due to the flat-fielding process 
(see details in the text, at the end of Section 2). The orientation of the north and east axes are shown. 



well as for instrumental transmission variations in the field- 
of-view using a flat-field observed at twilight. Next, we cre- 
ated wavelength-sorted data cubes corrected for airmass de- 
pendence, guiding shifts, cosmic rays, dark, and gain. We 
then corrected for gh osts, a common arte fact in FP datcQ, 
using the method of lEpinat et al.1 (l2008h . Fourth, we ap- 
plied a spectral smoothing using a Manning function, and 
then sky subtracted the data. We then combined all data 
cubes over all the cycles and nights observed to create a to- 
tal integrated cube. A spatial smoothing of 5 x 5 pixels was 
applied to allow better identification of diffuse gas. Finally, 
we fitted the line profiles with a single gaussian component 
and then corrected for the filter dependence, with a mean 
temperature of 10^ C. We were then able to extract the final 
velocity, dispersion and Ha integrated maps. 

For each galaxy, we show an image, summed over the 40 
channels, of the raw data in Fig. [2] (before applying the 5x5 
pixel binning and gaussian fits to the spectra). This figure 
shows that there is a contaminating reflection in the top and 
right side of the data (see rectangles outlined in red). This 
reflection arises from light reflected off the detector, and 
depended on the light present in the dome, and was there- 
fore brighter when taking the flats at twilight than when 
observing the galaxies. The data cubes were flat- corrected 
by dividing them with the normalized flat-field image. The 
continuum of the data cubes, in the regions of reflection, 
can therefore be inferior to that of the rest of the cube, and 
is not representative of the true continuum. This parasitic 



^ Ghosts are caused by reflections of bright sources (like a bright 
star or nucleus), and are located symmetrically opposite with re- 
spect to the optical axis. They will be at the same wavelength 
as the primary image, and can thus be well identified, both in 
position and in velocity. 



reflection only affected areas of diffuse emission. Section 3 
shows in more detail how this reflection was accounted for 
in the analysis of the diffuse extended emission. 



3 VELOCITY FIELDS 

We applied the reduction steps outlined in Section 2, and 
for each emission line we obtained the radial velocity mea- 
surement, standard deviation and integrated flux by fitting 
a gaussian function. These respectively gave the full res- 
olution velocity fields, dispersion maps and flux- calibrated 
Ha integrated maps shown in Fig. [3] (NGC 247) and Fig. g] 
(NGC 300). We also obtained lower resolution velocity fields 
(see Fig. [3] and Fig. |4l middle-right panels), by applying an 
additional gaussian spatial smoothing of wx = 55^^ where 
wx is the width of the smoothing function. These low resolu- 
tion maps were used to extract the kinematical parameters 
of the galaxies, which would have been difficult to do with 
the original full resolution velocity fields because of pertur- 
bations in the field (e. g. see velocity field of interarm region 
in NGC 247, where we were able to detect some gas). By 
using these less disturbed kinematical parameters, we then 
used the original full resolution velocity fields to extract the 
rotation curves. 

Excessive smoothing could induce a bias i n the determi - 
nation of the inclination towards lower values. ISosmal (|l978h 
determined that a rotation curve starts to properly trace the 
kinematics of a galaxy when the scale of the galaxy is more 
than 7 times the resolution of the data. Since our smooth- 
ing corresponds to less than 1/10 of the scale of the galaxies, 
the values we determine for the inclination should not be af- 
fected significantly. 
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Figure 3. Deep Fabry-Perot Ho; kinematical maps of NGC 247 obtained after a 5 x 5 pixels square binning and a single spectral 
gaussian fitting procedure, unless otherwise specified. Top-left: Ho; velocity field. Top-right: Dispersion map. Middle-left: Flux calibrated 
integrated map (scaled in 1/6 Rayleigh or lO"-*^^ ergs cm~^ s"-*^ arcsec"^). Middle-right: Ha velocity field obtained with an additional 
gaussian spatial smoothing of w\ = 55^^ Bottom: Position- velocity (PV) diagram along the major axis. The red and black lines are the 
rotation curves, taken along the major axis (with a slit width of 13 pixels) of a model velocity field. The black line is the rotation curve 
we derived for both sides, and therefore affected by the non-circular motions of the approaching side. The red line is the rotation curve 
we derived for the receding side alone, i. e. not affected by the non-circular motions. The latter fits much better the maximum intensities 
of the PV diagram. 



3.1 Annular ring binning for NGC 247 



Both 



Bland-H awthorn et alJ (|l997l ) and 

iHlavacek-Larr ondo et alJ (|201Q| ) succeeded in detecting 
diffuse ionized gas in NGC 253 by binning the data into 
annular rings. Since NGC 247 is also highly inclined 
{i ^ 74"^), we can apply the same kind of binning. The 



annular rings were centred along the major axis in the plane 
of the image and about the galaxy's dynamical centre. For 
a given opening angle, different regions corresponding to 
different radial intervals were binned together. We chose an 
opening angle of ^ 20° (i. e. 10^ on each side of the major 
axis), so that we would not overlap on too large a portion 
of sky, and made the binning intervals grow with radius. 
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Figure 4. Ha kinematical maps obtained for NGC 300, same caption as Fig. O Top-left: Ha velocity field. Top-right: Dispersion map. 
Bottom-left: Flux calibrated integrated map (scaled in 1/6 Rayleigh or lO"-*^^ ergs cm~^ s"-*^ arcsec"^). Bottom-right: Ha velocity field 
with a gaussian spatial smoothing of w\ = 55^^ Bottom: PV diagram along the major axis. The rotation curve shown in red is the one 
derived for both sides, and is the curve we adopted for the mass models. 



since diffuse emission becomes more difficult to detect as 
the radius increases (see Fig. [5]). 

As seen from our raw data in Fig. [21 the extended part 
of the northern arm of NGC 247 (beyond a radius of 11^) fell 
directly into the areas affected by the contaminating reflec- 
tions. Here, it became crucial to determine if any emission 
line detected was the result of diffuse emission or of the re- 
flection. The continuum observed in the affected regions was 
also not representative of its true value, and could amount 
to a negative value since it had been previously overcor- 
rected for during sky subtraction. However, the spectrum 
of the reflection did not vary significantly. By applying an- 
other annular ring binning, this time with opening angle of 
80^ to smooth sufficiently over the affected areas, we were 
able to obtain a typical response spectrum for the contam- 
inating reflection. If an emission line was only seen in the 
first (opening angle of 20"^) and not the second (opening 
angle of 80^) binning, then it was considered to be truly as- 



sociated with diffuse emission. Using this method, we were 
able to identify two extended radial intervals on the north- 
ern side of NGC 247 which had faint Ha emission (r = 11.2^, 
r = 13.1'). 



4 KINEMATICAL PARAMETERS AND 
ROTATION CURVE 

Four kinematical parameters must be determined from the 
velocity field to derive the rotation curve. They are the dy- 
namical centre (xo and ^o), the systemic velocity 14ys, the 
inclination i, and the position angle PA of the major axis. 
We use the gipsy program and ROTCUR task to derive them. 
These programs are based on a least square fittin g method 
that best reproduces the observed velocity field ([Begemanl 
Il989l ). 

Using the lower resolution velocity fields {wx = 55^'), we 
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Figure 5. Left: Velocity field of NGC 247 obtained after a 5 pixels x 5 pixels square binning procedure, annular ring binning, and a 
single gaussian adjustment to the emission profiles. The last two binned intervals (black contours) of the northern side show very faint 
and diffuse Ho; emission. Middle and right: Spectra of the last two binned intervals of the annular rings in order of radius. The star-like 
points are for the spectra of the diffuse emission (opening angle of ~ 20^) and the continuous lines are for the contaminating refiection 
spectra (opening angle of ~ 80^). The red vertical lines represent the location of the Ha emission. For the graph centred at a radius of 
13.1^, the parasitic refiection spectra (continuous line) was lowered in order to allow better comparison with the diffuse spectra (star-like 
points) by a value of 0.5 in units of 10"-*^^ ergs cm~-^ s"-*^ arcsec"^ A~^. As mentioned in Section 2 and Section 3, the continuum of the 
spectra is not representative of their true value. 



first find the dynamical centre and 14ys by keeping i and PA 
fixed. The initial values are chosen as those derived in the 
H I study (Carignan & Puche 1990 for NGC 247; Puche et al. 
1990 for NGC 300). To minimize errors due to deprojection 
effects, we exclude data in an opening angle of 60^ around 
the minor axis for NGC 247 and 35^ for NGC 300. We use 
a cosine-square weighting function on the remaining data of 
NGC 247, where each point is weighted by cos(^)^ with be- 
ing the angle from the major axis. This maximizes the weight 
of the points around the major axis. For NGC 300, since the 
inclination is smaller, we use a simple cos(^) weighting func- 
tion. 

We find V^ys = 164.8 ± 5.3 km s"^ for NGC 247 and 
of Vsys = 151.6 dz 4.2 km s"^ for NGC 300, in agreement 
with the systemic velocities found in the Hi studies (^-^161 
km s"^ for NGC 247, and 144.5±3.0 km s"^ for NGC 300). 
The dynamical centre we derive for NGC 247 coincides with 
the photometric one, but the one for NGC 300 was slightly 
shifted by ~ 30^^ {^^ 260 pc) towards the north-east. In the 
World Coordinate System, we find that the dynamical centre 
is located at a=00:54:50.4 ^=-37:39:48.0 for NGC 300, and 
at a=00:47:08.6 ^=-20:45:27.3 for NGC 247. 

Next, we determine i and PA, while the dynamical cen- 
tre and V^sys are kept fixed at the extracted values. Separate 
solutions are obtained for the receding and approaching side 
(i. e. using data points only on the receding or approaching 
side of the galaxy). For NGC 247, we find on average that 
i = 74.0 ± 5.8^ and PA= 171.0 ± 2.0^ and for NGC 300, 
i = 55.1 ± 8.0° and PA= -70.1 ± 2.9°. ROTCUR could not 
converge for a small number of radii, and the value of the 
parameters at these radii were interpolated using adjacent 
values. For NGC 300, the velocity field extends up to 13. 5^ 
but we were only able to extract parameters up to 12.0^ 
(even in the low resolution map, it was difficult to extract 
reliable kinematical parameters beyond 12.0^). To derive the 
complete rotation curve up to a radius of 13. 5^ we extrap- 
olated i and PA as the ones of the last data point beyond 



12. 0^ Some large fiuctuations were seen in the derived val- 
ues of i and PA, and to reduce their impact when deriving 
the rotation curve, we applied a median filter to the param- 
eters (median between 5 points). Fig. [6] and Fig. [7] shows 
the resulting values. Separate solutions are found for the re- 
ceding, approaching and combined sides. For NGC 300, the 
variations of i and PA at r > 9^ follow the same trend as in 
the H I data and are cons istent with the warp mentioned in 
ICarignan Puchd (Il990h . For NGC 247, more gas is seen 
on the northern (or approaching) side, and beyond a radius 
of 8.5^ only emission on the approaching side was found. 

Finally, using the full resolution velocity fields, we de- 
rive the rotation curves. Since NGC 247 does not show any 
sign of strong perturbations such as a warp or a bar, we de- 
rive the rotation curve using the average value of i and PA, 
i. e. we use the same value of i and PA throughout all radii, 
but derive the average value from the receding, approaching 
and combined sides separately. The average values of i and 
PA are shown in Fig. [6] with the coloured lines. To account 
for the warp in NGC 300, we must allow some variation of i 
and PA with radius, and therefore choose to derive the rota- 
tion curve using the median filtered values of i and PA, for 
each side separately. The bottom-left panels of Fig. [6] and 
Fig. [71 show the solutions found for the approaching and re- 
ceding side, and the bottom-right panels show the rotation 
curve derived for both sides, along with the proper error 
bars. The error bars are defined as the largest difference be- 
tween the rotation curve of the whole galaxy and that of one 
of the sides, or if larger, the intrinsic error determined by 
the ROTCUR task. The rotation velocities derived for both 
sides are given in Table [2] for NGC 247, and in Table g] for 
NGC 300. 

NGC 300 shows no signs of strong non-circular mo- 
tions. Thus, we simply use the rotation curve derived for 
both sides as the final adopted curve. However, NGC 247 
shows strong perturbations in the rotation curve between 
r ^ A' and r 8^ on the approaching side. We associate 
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Table 2. Derived Ho; rotation curve of NGC 247, when consid- 
ering both sides of the galaxy. Here, the large error bars do not 
represent the intrinsic uncertainty associated with the velocities. 
Instead, they reflect the large discrepancy arising between the so- 
lutions found for the rotation curve of the whole galaxy and that 
of the approaching side, the latter being affected by important 
non-circular motions. 
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these with non-circular motions of the interarm region on 
the approaching side of the galaxy (see Section 6.2.3 for a 
more detailed look at this region) . Applying a mass model to 
this distorted rotation curve proves to be very difficult, and 
does not give parameters representative of the true gravita- 
tional potential. Instead, we use the rotation curve derived 
for the receding side as the final adopted rotation curve of 
the galaxy, and include the two detections of faint Ha emis- 
sion seen in the extended radial intervals at r — 11.2' and 
r — 13. l^ Tableland the coloured points in the bottom- 
right panel of Fig. [6] show the values of the final adopted 
rotation curve. 

The velocities for the emission seen in the two extended 
radial intervals were derived using the kinematical parame- 
ters at maximum radius (r = 10.5^). For these points, along 
the X-axis, the error bar was chosen as the half width of the 
radial interval from which the emission line was extracted. 
Along the y-axis, we can not estimate the error bars in the 
same way as the other points (i. e. as the largest difference 
between the rotation curve of the whole galaxy and that of 
one of the sides or, if larger, the intrinsic error) since we 
only detect emission on the approaching side of the galaxy. 
The error bars were therefore chosen as the dispersion of the 
gaussian profile fitted to the emission line, since the largest 
uncertainty is on the location and shape of the profile, rather 
than on the kinematical parameters. An error on the order 
of 20^ for the inclination would be required to affect signifi- 
cantly the derived rotation velocities. Finally, the error bars 
associated with the velocities determined for the receding 
side (red points in the bottom-right panel of Fig. [6]), were 
taken as those derived by the ROTCUR task. 

In the bottom panels of Fig. [3] and Fig. |4l we show the 
position- velocity (PV) diagrams for each galaxy, along with 
the final adopted rotation curves in red, and for NGC 247, 
we include in black the rotation curve derived for both sides, 
which is affected by the non-circular motions. 



Table 3. Adopted rotation curve for NGC 247, including the 
detections of faint Ha emission associated with the radial intervals 
at r = 11.2' and r = 13.1^ (bold font). The error bars are taken 
as those derived by the ROTCUR task, except for those associated 
with the areas of extended emission (see Section 4). 
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5 MASS MODELS 

Observational data seem to favour fl at core mass models 
rathe r than cuspy core m o dels ( see Blais-Ouellette et al.l 



200ll: Ide Blok Bosm'3 l2002l: Ide Blok et al.1 l2003l : 
Kassin et al.1 l2006l : iSpano et al.l l20Q8l ). We therefore 



choose to model a dark halo in the form of an isothermal 
sphere, which has a density proflle as defined by Eq. 
(TJ This model is parameterized by a central density po 
and core radius Tc- The parameters are tied together by 
po = 9cr^/47rGrc^, where a is the one dimensional velocity 
dispersion. 



p{r) 



Po 



[1 + (r/rc)2]3/2 



(1) 



Based on a least-square fitting procedure, three param- 
eters must be determined during the mass model analysis 
(Carignan & Freeman 1985; Carignan 1985a). The first two 
are po and rc, which characterize the dark halo, and the 
third is the mass-to-light ratio (M/L)^, which characterizes 
the luminous component. We use the same s urface bright- 
ness pro file (B band) as inlPuche et al.1 (Il990l ) for NGC 300 
and . Carignan Puchel (Il990l ) for NGC 247 to describe the 
stellar disc. The absolute magnitude of the sun is set to 5.43 
in this band, and we take the same profile for the gaseous 
component as the one in the Hi papers, which includes a 
4/3 correction factor to account for helium. 

Best-fitting models generally reproduce well the rota- 
tion curves. However, various combinations of parameters 
can give the same results, i. e. these models are dege ner- 
ate (Ivan Albada Sancis il ll986l :[c arignan Fr eeman 19851 : 
button et al.ll2005h . Models such as those with no dark halo 
or those that are sub-maximum can be more interesting since 
they give only one solution. We applied a model with no dark 
halo, but this gave very poor results for both galaxies. 

Instead, we choos e to use the model developed by 
ICourteau Rixl (|l999l V These authors argued that for the 
majority of high surface brightness (HSB) late-type galax- 
ies, pure maximum stellar discs were not adequate. They 
then observed that the Tully-Fisher relation did not seem to 
show any dependence on surface brightness for HSB galaxies 
when plotted as Mr versus V2.2 (velocity at 2.2 exponen- 
tial scale lengths), and using this, they modelled the con- 
traction of dark haloes by adiabatic infall of baryons. Their 
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Figure 6. Top: Median filtered kinematical parameters of NGC 247 derived from the lower resolution velocity field for both sides (black 
symbols), i. e. solutions found considering data points along both receding and approaching sides, as well as separately for the receding 
(red symbols) and for the approaching (blue symbols) sides. The average values are shown with the horizontal lines. Bottom-left: Rotation 
curve derived using the average kinematical parameters of the top figures, with a systemic velocity of Vsys = 164.8 km s"-*^, and the full 
resolution velocity field. The continuous line, along with the derived error bars, represents the solutions found for both sides, while the 
receding and approaching sides are represented by the red and blue symbols, respectively. Bottom-right: Rotation curve derived for both 
sides, with the proper error bars (continuous line). These error bars reflect the large difference arising between the velocities derived for 
both sides, and those of the approaching side, the latter being strongly affected by non-circular motions. The y therefore do not represen t 
a large uncertainty in the intrinsic value derived for the velocity. The dashed line represents the Hi data from lCarignan Puchj (|l990l ). 
The coloured points illustrate the flnal adopted rotation curve, and include the two detections of diffuse Ha emission in blue. 



models then suggested that discs were sub- maximum with 
'^^disc/'^^totai 0.6, which requires that (M/L)^ be fixed at 
a specific value. For NGC 300, (M/L)^ must therefore be 
- I.8M0/L0 and for NGC 247, - 4.9Mq/Lq. For these 
particular galaxies, the ratios are consistent with the mod- 
els being maximum- discs. 

For NGC 300, we analysed a total of three rotation 
curves (see Fi^. [8] and Table [5]): first, the Hi curve from 
IPuche et alj (jlQQol l, second, the Ha curve derived here, and 
third, the combined Ha (r < 13.6^) and Hi (r > 13.6^) curve, 
since the latter extends further ou 10. For NGC 247 (see Fig. 
[9] and Table O , only the Ha and H i curves were analysed 
since our Ha rotation curve extends to the end of the Hi 
disc. For the H i cu rves, we include the e rror bars derived in 
IPuche et al.] (|l990l ) for NGC 300 and in lCarignan Puchd 
(I1990II for NGC 247, and for our Ha curves, we include the 



error bars associated with the adopted rotation curves (see 
Table 4 and Table 3, respectively). 



^ Since the i and PA are compatible between studies, we can 
combine both rotation curves. 



6 DISCUSSION 

We have obtained very deep Ha observations of NGC 247 
and NGC 300, in order to study the extent of the optical 
disc in these galaxies. For NGC 300, we detected emission 
all the way to the extremities of our field of view (24 x 24^), 
but suspect that there is still diffuse emission beyond. Other 
pointings and/or a larger field of view would be necessary to 
study the total extent of ionized gas in NGC 300. In the first 
section, we briefly discuss the results concerning this galaxy. 
We then address in greater detail the case of NGC 247, where 
we were able to detect Ha emission on the entirety of the 
H I disc. 
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Figure 7. Same as the previous figure (Fig. [6]), but for the galaxy NGC 300. The systemic velocity of Vsys = 151.6 km s"-*^ was used, 
and we used the median filtered values of i and PA to deriv e the rotation curv e, therefore allowing some variation with radius to account 
for the warp. The dashed line represents the Hi data from lPuche et all (|l990l V 



Table 4. Adopted Ho; rotation curve for NGC 300. The error 
bars are derived as the largest difference between the solutions 
found for the rotation curve of the whole galaxy and that of the 
approaching side or receding side, or if larger, the intrinsic error 
associated with the rotcur task. 
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NGC 300 is a late-type spiral gala xy with a ver y rich distri- 
bution of bright Hii regions ( Rizzi et al.ll2006l ). and a very 
extended st ellar disc reaching r — 2A' th at shows no sign of 
truncation (|Bland-Hawthorn et alJliooBI V The faintest level 
of the stellar disc reaches a surface brightness of 32 mag 



Table 5. Mass models for NGC 300 



Rotation Curve 


Hi^ 




36-cm ¥La 


H l+36-cm Ha 


Sub-Max 












1.8 




1.8 


1.8 


Dark halo^ 


Tc = 


173.6 


rc = 2.4 


rc = 9.7 




PO = 


4.9 


PO = 28.8 


PO = 7.4 



^ Puche et al. (1990) 

^ rc in kpc, a in km s"-*^ and po in 10~^ Mq pc~^ 

arcsec"^ (in V and B b and), and the stell ar population was 
found to be mostl y old (IVlaiic et al.[l2009l ). The Hi disc ex- 
tends to r ~ 20' (jPuche et al.lll99ol ). and seems to be flat 
beyond r = 15' (8 kpc). 

We have obtained very deep Ha observations of this 
galaxy. We reach EMs on the order of 0.1 cm~^ pc, and 
our rotation curve extends up to r = 13.5'. These are to 
our knowledge the deepest Ha o bservations of this ga laxy. 
Other Ha stu d ies in clude that of iMarcelin et al.l (|l985l ) and 
IHoodcs et al.l (|l996l ). The first obtained a rotation curve 
extending only to r = 10', and the second reached EMs 
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Figure 8. Maximum stellar disc mass models for NGC 300, where 
an isothermal sphere was used for the dark halo. The stellar disc 
is plotted with a small-dashed line, the gaseous component with 
a dot-dashed line, the halo with a long-dashed line and the total 
with the continuous line. Top: Mass model for the adopted H I 
rotation curve (diamond symbols) taken from lPuche et all (|l990[ ). 
Middle: Same but for the Ha rotation curve derived in this study 
(filled black circles). Hi data are also shown with diamond-shape 
points. Bottom: Same but for the combined Ha (r < 13.6^ or 
r < 7.1 kpc) + Hi (r > 13.6^ or r > 7.1 kpc) rotation curve. The 
parameters extracted are in Table O 



between 6 and 500 cm~^ pc. We detect significantly more 
diffuse emission. Fig. O (bottom-left panel) shows that the 
bright Hll refflons are emb edded in a large diffuse back- 
ground. iHoopes et aD (Il996l ) estimated that the diffuse frac- 
tion of Ha gas after scattering correction lies between 44 to 
54 per cent of the total Ha emission. Our study, which finds 
more diffuse gas, suggests that this fraction could be signifi- 
cantly larger. Since we were limited by our field of view, we 
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Figure 9. Same caption as Fig. [8] but for NGC 247. Here, we 
analysed the H I (top) and the Ha (bottom) rotation curve derived 
for the receding side (not affected by the non-circular motions), 
which includes the two detections of diffuse emission at large radii. 
The first is shown with diamond-like symbols, and the second is 
shown with the circles. The parameters extracted are in Table [G] 



Table 6. Mass models for NGC 247 



Rotation Curve 


Hl^ 


36-cm Ha (adopted curve) 


Sub-Max 






(M/Lb)^ 


4.9 


4.9 


Dark halo^ 


rc = 8.8 


rc = 5.0 




po = 8.6 


PO = 17.7 



« ICarignan Puchd (Il990h 

^ rc in kpc, a in km s~-^ and po in 10""^ Mq pc~^ 

suspect that there is more diffuse gas beyond our maximum 
radi us (r — 13.5^ ) . 

IPuche et al.l (|199Q[ ) found a distinct spiral structure on 
the southern-east side of the galaxy at r 20^ It would 
be very interesting to probe this structure, as well as to 
see if the Ha disc reaches the end of the H i disc, as it has 
been the case for other Sculptor group galaxies (Dicaire et 
al. 2008 for NGC 7793; Bland-Hawthorn et al. 1997 and 
Hlavacek-Larrondo et al. 2010 for NGC 253; and this study 
for NGC 247). 

In Fig. [71 we show the kinematical parameters and ro- 
tation curve we derive for NGC 300. Although our Ha data 
have a slightly better spatial resolution than the Hi data 
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(2.8'' X compared to ^" x ^" for the central and 17" x 12" 
for the outer parts in Hi), they could still be affected by 
beam smearing in the inner regions. Both rotation curves are 
nonetheless consistent, and the only significant difference is 
seen at large radii (r = 10 — 14') where the Ha rotation curve 
appears to drop and then raise again. This is seen in both 
the receding and approaching sides. The variation of the 
kinematical parameters with radius shows that there could 
be a warp in the disc starting at r ^ 9'. This is suggested 
by the continuous increase of the PA with radius sta rting at 
r ^ 9 , also seen in the Hi data (|Puche et al.lll99Ql ). mean- 
ing that the warp seems to affect both the Hi and optical 
disc. The PV diagram is shown in the bottom-panel of Fig. 
[H Here, we see that the rotation curve we derive seems to 
coincide with the maximum intensity at a given radius, in- 
dicating that the velocities we derive are able to reproduce 
the velocity structure in NGC 300. This diagram also con- 
firms that the dynamical centre is offset by 30" from the 
photometric one, which can be seen in the PV diagram as 
continuum emission at an angular offset of —30". 

We now briefly address the results concerning the mass 
models (see Fig. [8] and Table [5]). We applied a model which 
required (M/L)^ = LSMq/L©. This (M/L)^ is in good 
agr eement with the one pr edicted by stellar population mod- 
els. (Ben&^eJonS (|200lh predict for NGC 300, using a {B- 
V) magnitude of 0.59 mag given by the RC3 catalogue, that 
(M/L)^ in the B band is 1.2 Mq/L©. (M/L)^ = 1.8Mq/Lq 
is also consistent with the value required to maximise the 
contribution of the stellar disc in the inner regions (i. e. 
a maximum stellar disc model). However, this value of 
(M/L)^, and even the one predicted by stellar population 
models, seems to slightly overestimate the inner parts of the 
rotation curve, at least for two Ha velocity points. 

Fig. [8] shows that the maximum disc model reproduces 
well all three rotation curves, but the parameters derived for 
the Ha and combined Ha and H i rotation curves are signif- 
icantly different than those of the H i rotation curve alone. 
Our results suggest that the first two (Ha and Ha+Hi) 
have a smaller core radius, but this can be explained by the 
slightly lower velocities observed in Ha between r = 5 — 7 
kpc (r = 10— 14'). It seems that since the stellar component 
reproduces well the rotation curve at r < 3 kpc, the mass 
model then tries to adjust the dark halo component to repro- 
duce the shape of the curve for r > 3 kpc. Between r = 5 — 7 
kpc, the Ha rotation curve has lower velocities, which could 
make the model converge towards a more concentrated halo 
(smaller and more realistic rc). 



6.2 NGC 247 

6.2.1 Detection of extended Ha emission 

In NGC 247, we could study the extent of the Ha emis- 
sion on the entirety of the Hi disc. Our observations reach 
EMs of 0.1 cm~^ pc (not including the fainter emission 
detected using the annular ring binning in Section 3.1). 
This make s our study one of t he deepest Ha studies of 
NGC 247. iFerguson etaD (|l996l ) looked at diffuse ionized 
gas in NGC 247, but only reached a level of 6 cm~^ pc for 
Ha+[Nii], and a maximum radius of r = 10'. We detect 
much more diffuse emission, especially between the spiral 
arms (although still quite faint and with a disturbed veloc- 
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Figure 10. Flux calibrated integrated map of NGC 247 shown 
in the top panel. The bright Hll region to the north-west was 
believed to be an H ll ring causing a distortion in the rotation 
curve. We extract two regions along the north oppositely sy met- 
ric with respect to the major axis (shown with the blue and black 
contours), and plot in the bottom panel the derived projected 
rotation velocity as a function of radius. By projected rotation 
velocity, we mean that the observed radial velocity has been cor- 
rected for Vsys and z, but not for the cosine dependence on the 
azimuthal angle in the plane of the galaxy. The region within the 
blue contours contains most of the diffuse emission of the inter- 
arm region, and the region within the black contours contains the 
north H II region and some diffuse emission of the interarm re- 
gion. The bottom panel shows that the region dominated by the 
diffuse interarm emission has lower rotation velocities than when 
excluding this region. 

ity field i n the northern a rm) . We also find emission up to 
r = 13.1'. lFbrguson et al.l (|l996h determined that diffuse gas 
contributes to ^ 50 per cent of the total Ha luminosity in 
NGC 247. We suspect, as for NGC 300, that this fraction is 
much higher if we include all of the fainter gas we detect. 
This emphasises that diffuse ionized gas can contribute sig- 
nificantly in a galaxy, and that it can be used to probe its 

outer disc^ 

In iHlavacek-Larrondo et al ] (|2010h . we succeeded in 
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detecting [Nil] (A 6548 A) emission beyond the Hi disc 
in NGC 2 53, confirming the dec l ining part of tfie rota- 
tion curve. iBland-Hawtfiorn et alJ (| 19971 ) reported tfiat tfie 
[N ii] (A 6548 A) /Ha ratio in NGC 253 was quite strong (-1) 
beyond tti e H i disc, and so our detection of [Nil] was not 
surprising. iKennicutt eTaD (|2008l ) find a [N ii] (A 6548 A and 
A 6583 A)/Ha ratio of 0.48 for NGC 253, but only of 0.24 
for NGC 247, suggesting an even smaller fraction for [Nil] (A 
6548 A) /Ha and making it more difficult to detect this line 
in NGC 247. 

By binning the data in annular rings, we were able to 
detect extended faint Ha emission on the northern side of 
NGC 247 at r = 11.2' and r = 13.1'. These two detections 
have respectively an EM of ^ 0.02 cm~^ pc and ~ 0.05 
cm~^ pc. We estimate that these values are within a factor 
of two accurate, given that the mean absolute deviation of 
the difference between the annular ring spectrum and the 
contaminating reflection spectrum has a noise level of about 
40 per cent. Based on this estimate of the noise level, our 
detection at r — 11.2' represents a 3cr detection and the one 
at r = 13.1' represents at least a 2a detection. We also only 
detect Ha emission up to r = 8.5' on the southern side of 
the galaxy. Although it does seem that the northern side of 
NGC 247 has a l arger amount of gas comp ared to the south- 
ern side (see e. g. lCarignan Puchelll990l ). our observations 
could be biased because the northern side of the galaxy has 
typical radial velocities of 80 km s~^, corresponding to a 
filter transmission of 80 per cent (maximum value for our 
filter) , whereas the southern side has radial velocities of 250 
km s~^, corresponding to a filter transmission probability of 
only 65 per cent. However, one of the main conclusions of 
our observations is that, unlike NGC 253, we do not detect 
any optical emission beyond the Hi disc. We now examine 
the implications of this non detection in NGC 247. 

6.2.2 Implications of the extended Ha emission 

iBland -Hawthorn et al ] (|l997l ) proposed several explanations 
as to why NGC 253 had an optical disc larger than the 
H I disc. Among them, they first looked at the metagalactic 
ionizing background, which predicted an upper limit EM of 
0.25 cm~^ pc, but their Ha observations had EMs much 
higher. They then proposed photoionization models (Ly- 
man continuum, Lyc, ionizing photons), which predicted a 
much lower fraction of [Nil] /Ha beyond the Hi disc than 
their data suggested. Finally, they proposed ionizing pho- 
tons from hot young stars located near the central regions 
of the galaxy, which could look out towards the outer disc 
through the galaxy's warp. Through dust scattering, these 
photons would have to travel through less dense environ- 
ments, i.e. they would be less absorbed, and could ionize 
the gas beyon d the H i disc. Our [N ll] and H a observations 
of NGC 253 (iHlavacek-Larrondo et al]l2010l ) favoured this 
explanation. 

If this theory holds, this would mean that galaxies with- 
out a warp could not ionize gas beyond the H i disc. NGC 247 
is not warped, and we do not detect Ha emission beyond the 
H I disc. Hence, our results are consistent with this scenario. 

Our Ha measurements for NGC 247 fall below the pre- 
dicted upper limit of the metagalactic ionizing background. 
However, this theory predicts that there should be optical 
emission beyond the Hi disc, which we do not observe. As 



for the Lyc ionizing photons, UV GALEX maps of NGC 247 
show diffuse emission extending up to ^ 11-12', which co- 
incides with the first of the two annular rings where diffuse 
gas was detected. Since there is no significant UV emission 
beyond r = 13', this could also explain why we did not see 
optical emission beyond the H i disc. However, photoioniza- 
tion models still have many problems relating to the details 
of how these photons are able to escape to large distances. 
The exact ionizing source of diffuse ionized gas still remains 
unclear (see a recent review on the topic by Haffner et al.l 
l2009l ). but our non detection beyond the H i disc in NGC 247 
and detection beyond the Hi disc for the warped galaxy 
NGC 253 seem to favour the scenario in which photons are 
only able to leak out towards the outer disc through a warp. 

6.2.3 Kinematical analysis 

We show in Fig. O the velocity field, dispersion map and 
integrated map of NGC 247. Fig. 6 then shows that the 
rotation curve derived for the approaching (northern) side 
has a large drop in velocitie s at r ^ 5'. This feature was 
also seen in the Hi study by Carignan Puch e (1990), as 
well as the optical study by ICarignan Freeman (19851 ) , 
and was interpreted as being caused by non-circular motions 
associated with a bright Hii ring to the north-west. 

In Fig. 10, we explore in more detail this feature. The 
bright H ll region to the north-west at r 5' is believed to 
be the ring causing the shift in velocities. We extract two 
regions on the approaching side of the galaxy, oppositely 
symmetric with respect to the major axis. These are shown 
with the blue and black contours. The region within the 
blue contours contains most of the diffuse emission of the 
interarm region, and the region within the black contours 
contains the Hll region and some diffuse emission of the 
interarm region. In the bottom panel, we plot the derived 
projected rotation velocity as a function of radius for each 
region, which is the observed radial velocity corrected for 
the systemic velocity and inclination of the galaxy, but not 
for the cosine dependence on the azimuthal angle in the 
plane of the galaxy. Since both regions are symmetric with 
respect to the major axis, they will both be affected in the 
same way by the azimuthal angle dependence. The derived 
projected velocities can therefore be directly compared with 
one another. 

The bottom panel of Fig. 10 shows that the region dom- 
inated by the diffuse interarm emission (in blue) has much 
lower velocities. Hence, the Hii region to the north is not 
the cause of the drop in the rotation curve , which is not 
in agr eement with the interpretation of Cari gnan Puch3 
(|l990l V Instead, it seems that the diffuse region between the 
spiral arms to the north is causing the lower rotation ve- 
locities seen between r = 4 — 7'. The velocity field of the 
interarm region has a distorted pattern with radial veloci- 
ties varying between 50 — 250 km s""*^, suggesting that this 
region contains important non-circular motions. This is also 
supported by the PV diagram in Fig. [3l which shows diffuse 
emission at r ^ 4' (or r ^ 240 — 250") of large velocities 
that do not coincide with the maximum intensities. 

As mentioned earlier, in the PV diagram of Fig. [3l we 
show in black the rotation curve we derive for both sides, 
which is affected by the non-circular motions of the interarm 
region. This curve falls below the maximum intensities on 
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the receding side of the galaxy, consistent with the idea that 
it does not represent the true rotation within the galaxy. In 
red, we show the rotation curve derived only using the reced- 
ing side, which follows well the maximum intensities in the 
PV diagram, and therefore represents well the true rotation. 
The latter was used when deriving the mass models. 

Fig. [9] and Table [6] show the solutions derived for the 
mass model analysis. The model required a (M/L)^ ratio of 
4.9 M0/L0, and gave a good fit for both rotation curves. 
The (M/L)^ ratio is qui te high when compar ed to stellar 
population models (e. g. iBell de JondboOll ). which pre- 
dict a ratio of 1.0 M0/L0 in the B band, using a (B-V) 
magnitude of 0.56 mag as given by the RC3 catalogue. How- 
ever, it is interesting to note that (M/L)^ = 4.9 M0/L0 is 
again consistent with the value required to maximise the 
contribution of the stellar disc in the inner regions, just like 
in NGC 300. Our results show that the Ha rotation curve 
requires a larger po (halo is more concentrated). This could 
be because of the slightly higher Ha velocity points between 
r = 3 — 7 kpc, and the two extended points found with the 
annular binning method which have lower velocities than 
the H I rotation curve. The model could therefore be trying 
to find a more concentrated dark halo. 



7 CONCLUDING REMARKS 

This paper is the third in a series that aims to study dif- 
fuse ionized gas of the Sculptor group galaxies using very 
deep optical observ ations. The first was carried out by 
iDicaire et all (|2008l ) and confirmed that NGC 7793 had 
a truly de clining rotation curve. The s econd looked at 
NGC 253 (iHlavacek-Larrondo et alll2010l ) and found very 
extended [N li] emission beyond the H I disc, showing a sig- 
nificant decline in the rotation curve. The third is presented 
here, and looks at Ha emission in NGC 247 and NGC 300. 
For NGC 247, we have detected diffuse Ha on the entirety 
of the Hi disc, but not beyond. For NGC 300, Ha emission 
is detected on the entirety of the allowed spatial coverage, 
but since we are limited by our field of view, we suspect that 
there is more optical emission beyond. However, our study 
still remains the deepest Ha study of these two galaxies, 
and from there, we have obtained the most extended opti- 
cal rotation curves thus far for both galaxies. To account 
for the existence of diffuse optical emission beyond the Hi 
disc in the warped galaxy NGC 253, but not in NGC 247 
(which is not warped), our observations favour ionization 
by hot young stars in the central region of the galaxy, which 
only through a warp, can let photons escape and ionize the 
interstellar medium in the outer parts. 
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